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We report a high energy-resolution neutron backscattering study, combined with in-situ diffrac- 
tion, to investigate slow molecular motions on nanosecond time scales in the fluid phase of phospho- 
lipid bilayers of l,2-dimyristoyl-sn-glycero-3-phoshatidylcholine (DMPC) and DMPC/40% choles- 
terol (wt/wt). A cooperative structural relaxation process was observed. From the in- plane scatter- 
ing vector dependence of the relaxation rates in hydrogenated and deuterated samples, combined 
with results from a 0.1 /is long all atom molecular dynamics simulation, it is concluded that cor- 
related dynamics in lipid membranes occurs over several lipid distances, spanning a time interval 
from pico- to nanoseconds. 

PACS numbers: 87.16.dj, 87.14.Cc, 83.85.Hf, 87.15.ap, 83.10.Mj 



It is speculated that atomic and molecular motions in 
regions of biomolecular systems with strong local inter- 
actions are highly correlated on certain range of time 
and length scales [H [2]. In proteins intra-protein cor- 
relations are believed to be essential for their biological 
functioning, such as protein folding, domain motion and 
conformational changes. Very recently, inter-protein cor- 
relations in protein crystals and also membranes have 
been reported from experiment and simulation [3J |J . Ex- 
perimental and computational effort has been invested 
to study collective molecular motions in phospholipid 
model membranes [3 [6l [7l |8] to understand the possi- 
ble impact on physiological and biological functions of 
the bilayers, such as transport processes [5], and eventu- 
ally their implication on function of membrane-embedded 
proteins. While fast (picosecond) propagating collective 
microscopic fluctuations in the plane of the bilayer can be 
understood as sound waves [TO] IH], the slow (nanosec- 
ond) in-plane mesoscopic fluctuations (undulations) are 
governed by the elasticity properties of the bilayers [T^ . 

We studied dynamical modes at nearest neighbor dis- 
tances of the lipid molecules using the neutron backscat- 
tering technique [T3]. These modes are too fast to be 
accessed by x-ray photon correlation spectroscopy and 
the lateral length scales involved are too small to be re- 
solved by dynamic light scattering or the neutron spin- 
echo technique. Selective deuteration was used to dis- 
criminate relaxations due to collective molecular motions 
from relaxations arising from localized, single molecule 
excitations. In this Letter, we examine results of inelas- 
tic neutron scattering experiments that demonstrate the 
existence of slow coherent motion of lipid molecules in 
the fluid phase of phospholipid bilayers. From the in- 
plane scattering vector dependence (g||) of the measured 
relaxation rates, combined with results of a 0.1 /is long all 
atom molecular dynamics (MD) simulation, we flnd that 
the cooperative structural dynamics in lipid membranes 
occurs over several lipid distances, spanning a time inter- 



val from pico- to nanoseconds. 

The experiments were carried out at the cold neu- 
tron backscattering spectrometer IN16 [14] at the Insti- 
tut Laue-Langevin (ILL) with an energy resolution of 
about 0.9 /leV FWHM (A=6.27 A). An energy trans- 
fer of —15 ^eV< E < +15 /ieV and a q range of 
0.43 A-i < q <1.92 A-i were scanned, accessing 
time scales of 0.28 ns<t<4.6 ns and length scales of 
3.2 A<d<14.6 A. A separate line of diffraction detectors 
allowed to determine the bilayer structure in-situ. Pro- 
tonated (DMPC-h) and partially (acyl chain) deuterated 
l,2-dimyristoyl-sn-glycero-3-phoshatidylcholine (DMPC- 
d54), along with protonated cholesterol were obtained 
from Avanti Polar Lipids. Highly oriented multi lamellar 
membrane stacks of several thousands of lipid bilayers of 
DMPC and DMPC/cholesterol were prepared by spread- 
ing a solution of typically 25 mg/ml lipid and lipid/40% 
cholesterol (wt/wt) in trifluoroethylene/chloroform (1:1) 
on 2" silicon wafers, followed by subsequent drying in 
vacuum and hydration from D2O or H2O vapor. Twenty 
such wafers (flve for DMPC-h) separated by small air 
gaps were combined and aligned with respect to each 
other to create a "sandwich sample" consisting of several 
thousands of highly oriented lipid bilayers (total mosaic- 
ity about 0.5°). The deuterated (protonated) samples 
had a total mass of about 400 mg (100 mg) of lipid or 
lipid/cholesterol. The samples were mounted in a her- 
metically sealed aluminum container within a cryostat 
and hydrated from D2O or H2O vapor. By aligning the 
bilayer normal at 135, respective 45 degrees with respect 
to the incoming neutron beam, the momentum transfer 
could be placed in the plane of the bilayers (q^^) or per- 
pendicular to the membranes (qz). The lamellar spac- 
ing of the DMPC sample was determined to (iz=54 A 
at T=30°C, which corresponds to a relative humidity of 
RII=99.6% p[5]. Five different samples (referred to as 
SI,..., S5) were prepared: 
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Neutrons are scattered by the atomic nuclei, and each 
element intrinsically has non zero coherent and incoher- 
ent scattering cross sections. In experiments the differ- 
ent contributions can be enhanced with respect to each 
other, by, e.g., using different isotopes or selective deuter- 
ation. While in protonated samples the incoherent scat- 
tering is dominant and the time-autocorrelation function 
of individual scatterers is accessed, (partial) deuteration 
emphasizes the coherent scattering and probes the pair 
correlation function. By careful analysis of the scatter- 
ing of all samples SI to S5, it was possible to emphasize 
between the collective and self motions in the underlying 
dynamics. Sample S5 allowed us to investigate the effect 
of cholesterol incorporation on the collective dynamics 
of the lipid backbone because protonated cholesterol was 
used. 

Figure [ija) depicts in-plane diffraction patterns, i.e., 
structure factors S{q\\), at 303 K. Subtracting S{q\\) for 
samples S3 and SI emphasizes the coherent scattering of 
the lipid acyl chains. Similarly, subtraction of S{q\\) for 
S3 and S5 emphasizes the coherent lipid tail scattering 
in the cholesterol sample. Finally, to emphasize coherent 
water scattering, the diffraction pattern of S4 was sub- 
tracted from that of S3. There are two peaks in >5'((7||) 
of the lipid chains around g|| = 0.7 and 1.4 A^^, 
which correspond to the nearest neighbor distances of 
phospholipid head groups and the acyl chains in DMPC, 
respectively. The peaks shift to smaller g|| values in the 
presence of cholesterol, indicating a slight increase in the 
nearest neighbor distances. The water scattering is en- 
hanced towards the water correlation peak around 2 A~^, 
but there is no pronounced diffraction peak within the in- 
strumental resolution. 

Figure [2] exemplarily depicts spectra measured at g|| = 
0.5 A^^ and 1.5 A^^ for SI and S4 at a temperature 
of T=303 K, in the fluid phase of the bilayers. A fast 
and a slow relaxation process were visible in the data 
and data were therefore fitted using two Lorentzians and 
convoluted with the instrumental resolution. Because 
the faster process (the broad Lorentzian) could not be 
determined with sufficient accuracy, this Letter focuses 
on the slower relaxation process (^^-"1). Wi might con- 
tain contributions from incoherent and coherent scatter- 
ing. A quantitative determination of the coherent and 
incoherent relaxation rates would require purely coher- 
ent and purely incoherent scatterers or neutron polar- 
ization analysis (which has not yet been realized on a 
neutron backscattering spectrometer) to discriminate the 
different contributions to wi. A convolution of incoher- 
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FIG. 1: (Color online), (a) Diffraction patterns generated 
by subtraction of S{q\\) for pairs of samples S1,...,S5, which 
emphasize contributions from collective molecular motions of 
the lipid acyl chains, and the water molecules (see text for ex- 
planation), (b) Aoji for the lipid acyl chains, the membrane 
hydration water, and the lipid chains 'distorted' by incorpo- 
rating 40% cholesterol. For the water, the difference of Acji 
of samples SI and S2 was taken. 



ent and coherent relaxation processes, each described 
by a Lorentzian peak shape (implicating a single ex- 
ponential relaxation process), was assumed. Since the 
convolution of two Lorentzians with widths Auj{ and 
Aw" is again a Lorentzian peak shape with a HWHM of 
Awi = Au'i+Alui, by subtracting the Lorentzian widths 
of spectra of protonated (mainly incoherent scatterers) 
from the (partially) deuterated samples (mainly coher- 
ent scatterers) , it was qualitatively possible to distinguish 
the cooperative dynamics relative from the single-particle 
dynamics. 

Fig. [ijb) shows the difference in the quasielastic width 
Aci;i(g||) for the slow process for the hpid acyl chains, 
the membrane hydration water, and the lipid chains 'dis- 
torted' by incorporating 40% cholesterol. So does, e.g., 
subtracting the widths of spectra S1-S3 emphasize the 
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FIG. 2: (Color online). Exemplary inelastic spectra in the 
fluid phase of the bilayers. Left column: = 0.5 A~^; right 
column: g|| = 1.5 A-\ a) DMPC-d54/D20 (Sample SI), b) 
DMPC-h/H20 (Sample S4). Data have been fitted using the 
instrumental resolution (dotted) , two Lorentzian peak shapes 
(dashed), and a constant background (dashed). The total fit 
is shown as solid line. 
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FIG. 3: The slow relaxation times Tc (closed circles) and Ts 
(open squares) as a function of wave- vector q\ | . Also shown 
is the static structure factor 5(51 1) (open triangles). Inset: 
Is{q,t) (open symbols) and Ic{<l,t) (filled symbols) for g|| = 
2.53 (lower curves) and 1.42 (upper curves). 



coherent hpid tail dynamics, analogous to the diffraction 
discussed above. Awi((7||) shows a pronounced minimum 
around q^^ —1.22 for the lipids and lipids/cholesterol 
sample and qualitatively reproduces the behavior of 
S{q\\) in Fig. [T|a). This points to a strong coupling of 
dynamical properties to the static structure of the system 
and to a softening of the coherent relaxation rate around 
the maximum of the static structure factor S{q\\). The 
slowing down of the decay of the density autocorrelation 
function for distances corresponding to nearest neighbors 
is known as de Gennes narrowing |16j . The momentum 
dependence of the dynamical parameters in the collective 
dynamics case is not trivial as in the single-particle case 
but points to coherent structural relaxation. Because 
the out-of-plane component {qi_) of the scattering vec- 
tor (also in Fig.[l|b)) shows no minimum, it is concluded 
that the relaxation process is confined to the plane of the 
membranes. 

To gain further insight in the above relaxation process, 
we performed an all atom MD simulation of DMPC and 
analyzed the motion of the carbons within the lipid tails. 
The fully solvated system, built from a pre-equilibrated 
DMPC bilayer [T7], contained 128 lipid and 2577 TIPS 
water molecules. The MD simulation was performed 
with NAMD-2.6 [18] using the CHARMM27 [I9] force 
field. The pair interactions were turned off smoothly from 
10 A to 12 A, and the long-range electrostatic interac- 
tions were calculated using the PME method [5D]. The 
simulation was performed in the NVT ensemble using 
periodic boundary conditions. A constant temperature 
of T=303 K was maintained by employing a Langevin 
thermostat with a coupling constant of 0.05 ps""'^. After 
proper energy minimization and thermal equilibration, a 
production run of 0.1 fis was carried out on 40 CPUs 



of a dual core 2.8GHz Intel Xeon EM64T cluster with a 
performance of 0.2 day/ns. 

To examine the structural relaxation of the lipid tails 
we calculated the in-plane incoherent (self) Is{q\\,t) — 
(1/^) ^„(e"*?"i|-[*'"(°)-^"'"(*)l) and the coherent Ic{q\\,t) = 

(1/^) E„ X;m(e~''ii'^""^°^"*''"^*^') intermediate scattering 
functions. The sums are taken over the N carbon atoms 
in the lipid tails and r„(i) is the position of the n*'^ atom 
at time t. The scattering vector q^^ is parallel to the x — y 
plane of the lipid bilayer. By examining Is/c{q\\yt) we 
can isolate purely coherent and incoherent contributions 
to the scattering due to the carbons within the lipid tails. 

The minimum in AaJi(g||) shown in Fig. |l] suggests a 
longer relaxation time for Ic{q\\,t) than for Is{q\\,t) for 
qy values around the peak of S{q\\) = lc{q\\,0). Shown 
in inset to Fig. [s] are Is/c{q\\,t) for q\\ = 2.5 A~^ and 
1.42 A"^ For gy = 2.5 A~^ the decay times for Is{q\\,t) 
and /c((7j[,t) are almost identical, but /c((7||,t) decays 
slower for g|| around the first peak of the static struc- 
ture factor. We defined the decay time r((7y) as when 
Is/c{q\\,T{q\\))/Is/c{q,0) = Shown in Fig. [s] are 

Ts{q\\) and Tc(gy) along with 5'(gy). Note that the non- 
monotonic behavior of Tc{q\\) yields a minimum in the 
difference Tg — around the peak of >S'(gy). Thus, the 
MD simulation results are in qualitative agreement with 
those from the experiment and reinforces the conclusion 
that the difference of the quasielastic widths Awi(q||) is 
due to correlated dynamics of the lipid acyl chains. 

To determine the temporal and spatial extent of the 
correlated motion of the lipid tails, we calculated the 
correlation function 
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FIG. 4: Measures of the correlated dynamics T{t) = 
J G{r, t)dr. r{t) is a measure of the size of the correlated 
displacement at a time t, and G{r, t) measures the spatial 
range of the correlations. 
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where Ar„(<:) = |^ri,||(i) — ^n.||(0)| represents the displace- 
ment of caxbon atoms in the x ~ y plane of the mem- 
brane, and the prime denotes a restricted double sum in 
which carbon pairs belonging to the same lipid tail are 
excluded. In the absence of displacement correlations, 
ffA(r,t) = .91 (r) = {VlN^){Y.r.,rr!5{T-\rr.-r^\)). Thus 
= /[saI?", i) /51(f) — 1] dr provides a measure of the 
correlation of in plane displacements of the carbon atoms 
within the lipid tails between different lipids [21 . 

Examination of r(t), Fig.Qa), shows that the displace- 
ments of the lipid acyl tails are correlated for times from 
around 1 ps to 1 ns, i.e., for much longer than the micro- 
scopic collision time but shorter than the time needed for 
a lipid to diffuse one lipid diameter. This is exactly the 
time window of the structural relaxation probed by the 
neutron scattering experiments and corresponds to the 
decay times t^/^. Thus, we associate Acji observed in 
the experiments and the difference in decay times Tc — 
to correlated motion of lipid tails. 

Finally, the spatial decay of G(r, t) = {r, t) / gi (r) — 1 
yields the length scale associated with the correlated dis- 
placements. Shown in Fig.jlja is G{r, t) for ti — 2.5 ps (at 
the peak of r(t)) and t2 = 1-1 ns. For both times G(r, t) 
does not decay to zero until around 30 A. Therefore, the 
lipid tails displacements are correlated for at least four 
to six lipid diameters between ti and ^2- 

In conclusion we found experimental evidence for a co- 
operative structural relaxation process in fluid phospho- 
lipid membranes. An all atom molecular dynamics simu- 
lations demonstrated that the displacements of the lipid 
tails are correlated for up to six lipid diameters for a time 
span from picoseconds to around a nanosecond. Experi- 
mental data of the sample containing cholesterol proved 
that the process is present not only in pure lipid bilay- 
ers, but also in more complex and relevant membranes. 
A possible implication of this motional coherence of the 
lipid acyl chains is, e.g., that information about a lo- 



cal structural perturbance of the lipids could propagate 
in the bilayer, which might be relevant for the under- 
standing of processes and functions involving collective 
structural changes. 
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